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Abstract 

The electromagnetic form factors of the nucleon are calculated at order q 4 in the relativistic 
formulation of baryon chiral perturbation theory. In order to obtain a consistent power counting 
for the renormalized diagrams we make use of the extended on-mass-shell renormalization scheme 
which we have discussed in a recent paper. We analyze the Dirac and Pauli form factors, F\ and 
i*2i as well as the electric and magnetic Sachs form factors, Ge and Gm- Our results are compared 
with those obtained in the heavy-baryon approach and in the (relativistic) infrared regularization. 
In addition, the Sachs form factors are compared with experimental data. We confirm previous 
findings that a one-loop calculation of the electromagnetic form factors does not generate sufficient 
curvature for the Sachs form factors G P E , G P M , and G^. Moreover, the electric form factor of the 
neutron is very sensitive to higher-order contributions. 
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I. INTRODUCTION 



Chiral perturbation theory (xPT) for mesons 0, 0| is the effective field theory of low- 
energy QCD in the vacuum sector (baryon number zero), expressed in terms of the effective 
degrees of freedom active at low energies. These are the Goldstone bosons resulting from the 
spontaneous symmetry breaking of the chiral group SU(A^)^ x SU(iV)#, reflecting the sym- 
metry of the QCD Lagrangian for N massless quark flavors, down to its vectorial subgroup 
SU(iV)y, denoting the symmetry of the vacuum. The effective Lagrangian is organized in 
terms of a simultaneous expansion in powers of (covariant) derivatives and quark-mass terms, 
and a perturbative expansion of observables is expected to work for momenta sufficiently 
small in comparison with the scale of chiral symmetry breaking, A x ~ 1 GeV. Perturbative 
calculations in the mesonic sector are feasible due to the existence of a consistent power 
counting scheme 0, 0| establishing a connection between the chiral expansion and the loop 
expansion (see, e.g., Ref. [j| for a recent overview). 

The one-nucleon sector (baryon number one) turns out to be more complicated ji| because 
of the presence of a new mass scale — the mass of the nucleon — which does not vanish in the 
chiral limit. In particular, loop diagrams containing internal nucleon lines may contribute to 
low orders in the chiral expansion. Therefore, the relation between the chiral expansion and 
the loop expansion appears to be lost 0]. This problem has been solved in the framework 
of heavy-baryon chiral perturbation theory (HB;\PT) jH[ 0], where the power counting is 
restored in terms of an additional 1/m expansion, and perturbative calculations are again 
possible. However, HB;\PT has its own shortcomings: the corresponding perturbation se- 
ries fails to converge in part of the low-energy region [jj. Moreover, at higher orders in the 
chiral expansion, the expressions due to 1/m corrections of the Lagrangian become increas- 
ingly complicated j8|,|9|. These disadvantages are related to the nonrelativistic momentum 
expansion in this approach. 

Recently, various methods have be en p roposed to establish a consistent power counting 
within a relativistic approach [HI 11 , E2 EH, EH, EH, US 



171. In this work, we will discuss 



the electromagnetic form factors of the nucleon within the extended on-mass-shell (EOMS) 
renormalization scheme of Ref. [ItJ • In this scheme renormalization is performed by expand- 
ing a given Feynman diagram in terms of small quantities and subtracting those terms which 
violate the power counting. Since the subtracted terms are regular in the small quantities, 
they can be absorbed in a finite number of low-energy constants of the chiral Lagrangian. 
So far, the new method has been applied to a discussion of the mass and the scalar form 
factor of the nucleon at 0(q A ) 

The electromagnetic form factors parameterize the single-nucleon matrix element of the 
electromagnetic current operator and thus reflect the electromagnetic structure of the nu- 
cleon. The matrix element is usually described in terms of either the Dirac and Pauli form 
factors, Fx and F 2 , or the electric and magnetic Sachs form factors, Ge and Gm- The latter 
are popular, because their Fourier transforms in the Breit frame come closest to an inter- 
pretation as the distribution of charge and magnetization inside the nucleon |l8L Eil | (for 
a recent discussion and an overview of the existing relevant form factor measurements, see 
Refs. jiol . 2ll|). For low momentum transfers, Q 2 < 1 GeV 2 , the Sachs form factors G P E , G\ 
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and G^ are reasonably well described by a dipole form factor Gd, 

G p e « « « Gd , 
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where Gd(Q 2 ) = [l+Q 2 / (0.71 GeV 2 )]~ 2 . As has been stressed in Ref. [21], there is essentially 
no physical motivation for such a form and it is also not able to describe "fine structure" 
in the data. The electric form factor of the neutron, G%, is the least precisely known 
due to the lack of a free neutron target and the fact that it is small. However, recent 
polarization experiments have improved our empirical knowledge of G% considerably (see 



Ref. j2l| and references therein). Clearly, a quantitative description of these four form 



factors is a stringent test for any theory or model of the strong interactions (see, e.g., Refs. 
H 0, H H ESJ). In chiral perturbation theory, these form factors have 
been calculated within the early relativistic approach p[, HB%PT |(| 0], the small-scale 
expansion j^, and within the relativistic infrared regularization approach |33j | . 

Our work is organized as follows. In Sec. |H] we shortly review the effective Lagrangian 
and list the interaction terms relevant for our calculation. In Sec. IIHI the Dirac and Pauli 
as well as the Sachs form factors are discussed. Our results are compared with the HB^PT 
calculation, the method of infrared regularization, and experimental data. A summary and 
some concluding remarks are given in Sec. II VI Some technical details can be found in the 
appendices. 



II. CHIRAL LAGRANGIAN 

The relevant effective Lagrangian contains two pieces, 

where the first part is the pure mesonic Lagrangian and the second part consists of terms 
which are bilinear in the isospin doublet $ containing the proton and neutron fields. The 
different terms of the effective Lagrangian can be organized according to the chiral derivative 
and quark- mass expansion |l|, |2j, |4| , 

where the subscripts (superscripts) in C w (C^n) refer to the order in the expansion, and the 
ellipsis denotes terms of higher order which are not relevant for our calculation. 

From the mesonic sector we only need the lowest-order Lagrangian [0(g 2 )] including the 
quark-mass term and the coupling to an external electromagnetic field A? in terms of the 
isovector vector field = -ef A* (e 2 /4n w 1/137, e > 0) 0, 1 

F 2 / \ F 2 M 2 / \ F 2 r i 

C 2 = — Tr (dfjUd^) + Tr [U^ + U)+ i— Tr [(d^UW + d^U^] + ■ ■ ■ . (2) 

Here, F denotes the pion decay constant in the chiral limit, F n = F[l + 0(rh)} = 92.4 
MeV, and M 2 = 2Brh is the lowest-order prediction for the pion mass squared, where B 
is related to the quark condensate (qq)o in the chiral limit 0, 35[. We make use of the 



isospin-symmetric limit m u = = rh. The triplet of pion fields is collected in the matrix 



1 In the mesonic sector, the isoscalar coupling to the external electromagnetic field only becomes relevant 
at 0(q 6 ) Q. 
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u. 




%/27T — 7T 

The lowest-order pion-nucleon Lagrangian is given by 



(3) 



with 



u. 



In Eq. (J3J), m and <^ refer to the chiral limit (at fixed strange-quark mass) of the physical 
nucleon mass and the axial- vector coupling constant, respectively. In the definition of the 
covariant derivative we follow Ref. Q, where T M only contains traceless external fields and 
the coupling to the isosinglet vector field is considered separately 

The next-to-leading-order pion-nucleon Lagrangian contains seven low-energy constants Cj 
H|, where only the terms proportional to C\, C2, C4, cq, and C7 contribute to the calculation 
of the electromagnetic form factors at 0(q 4 ), 



Cl Tr( X+ )M - 
.c 4 



c 2 

4m 2 
c 6 



*lr( V „)DW + h.c. 



(4) 



where h.c. refers to the Hermitian conjugate and 

X+ = M 2 (C/ + C/t), 



± 



,9 / _ -9 1 _ ? - 17 71 
- £^77, - % [r M ,r„] , 



with r M = Z M = — and # = — eA^/2. With the convention of choosing /+ to be 

traceless, C6 (07) will be related to the isovector (isoscalar) magnetic moment of the nucleon 
in the chiral limit. 

The complete Lagrangians at third and fourth order can be found in Refs. jH E| and |j| , 
respectively. Again, we only list the terms needed for the calculation of the electromagnetic 
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form factors, 2 



£ (3) 



-d^ 



D",f, 



D"* + h.c. 



2m 

+—d 7 ^(d^)D v ^ + h.c. + ■ ■ 



in 



-2e 54 - - e 74 [£>\ [D X , f+ }] 



-2 eios v$Tr (x+) - - e W6 /+ Tj (x 
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The values for the low-energy constants (LECs) have to be determined from empirical input. 
In general, their numerical values depend on the renormalization scheme in question. As we 
will see later, some of the constants can be fitted to the anomalous magnetic moments and 
to the charge and magnetic radii of the nucleon. 



III. ELECTROMAGNETIC FORM FACTORS OF THE NUCLEON 
A. Definition of the electromagnetic form factors 

In terms of QCD degrees of freedom the interaction with an external electromagnetic 
field Afj_ is given by 

where the electromagnetic current operator reads 

J"(z) = -u{x)^u{x) - ~d(x)-fd{x) + ■■■ = q(x) (- + ^) Yq{x) + •••■ (5) 

The ellipsis denotes the contribution due to the heavier quarks which we do not consider 
here. The electromagnetic form factors are defined via the matrix element 



(N(p f )\r(0)\N(p t ))=u(p f ) 



u(pi), N = p,n, (6) 



where q = Pf — Pi is the momentum transfer and Q 2 = —q 2 = — t > 0. 3 The functions 
{Q 2 ) and FJ?(Q 2 ) are the Dirac and Pauli form factors of the nucleon, respectively. At 
Q 2 = 0, these form factors are given by the electric charges and the anomalous magnetic 
moments in units of the charge and the nuclear magneton, respectively: 

Ff (0) = 1, (0) = 0, iff (0) = Kp = 1.793, F 2 n (0) = K n = -1.913. (7) 

In that sense, Eq. (JHJ) provides a "natural" extension of the electromagnetic vertex of a 
"point particle" with an anomalous magnetic moment. 4 In the actual calculation, it is more 



2 We took the Lagrangians of Ref. [ij and made the replacements F+ — > /+ and Tr(F — > Avff . 
Since we discuss the form factors in the space-like region, here we adopt the convention of taking Q 2 = —t 
as the argument of the form factors as is common practice in the context of electron scattering. 

4 For a discussion of the implications of gauge invariance in its strong form and alternative forms of the 
vertex, see, e.g., Ref. [361 . 
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convenient to work in the isospin basis 

F^ = F! + F^ Ft ] = F?-F?, i = l,2, (8) 
so that the electromagnetic form factors are obtained as follows, 



Fr = ^ s) + T jFt\ . = 1,2. 



Experimental data are usually analyzed in terms of the electric and magnetic Sachs form 
factors [3 [13, (Q 2 ) and G^(Q 2 ), defined by 



Q_ 2 

4m? N 



g e(Q 2 ) = F \{Q 2 ) ~ 7^tF 2 n (Q*), (9) 



G^(Q 2 ) = F^Q^ + F^Q*). (10) 

The Fourier transforms of the Sachs form factors in the Breit frame can be related to the 
distribution of charge and magnetization inside the nucleon. For a recent discussion of the 
limits of such an interpretation, see Ref. jiol ]. 



B. Calculation of the Dirac and Pauli form factors 

The graphs contributing to the electromagnetic current matrix element of Eq. (JHJ) up 
to and including 0(q 4 ) are shown in Fig. ^ (we do not display external leg corrections). 5 
In particular, the loop graphs labeled by the numbers (5) to (9) in Fig. ^ contribute at 
0(q 3 ) while the graphs (10) to (12) contribute at 0(q 4 ). The individual unrenormalized 
contributions to the form factors F\ and F2 can be found in App. El To renormalize the 
form factor diagrams we first apply the subtraction scheme used by Gasser jmd Leutwyler 
0, H| which we denote by modified minimal subtraction scheme of %PT (MS). We then 
perform the additional finite subtractions according to our EOMS scheme (see Ref. [l7| for 
details). For this we determine the chiral order of a given diagram by applying the standard 
power counting and then subtract those terms which violate the piower counting. Such finite 
subtractions are achieved by expanding the quantities of the MS scheme in terms of the 
parameters of our EOMS scheme. This expansion produces the counterterm contributions 
which are responsible for the required additional subtractions. In view of Eq. (jUJ), an 0(q A ) 
calculation of the electromagnetic form factors yields F\ to 0(q 3 ) and F 2 to 0(q 2 ), since 
both the polarization vector e M and the four-momentum q^ of the virtual photon count as 
0(q). Diagrams potentially violating the power counting are loop diagrams with internal 
nucleon lines. In the present case, we find that the diagrams (5), (8), and (10) of Fig. [T]need 
subtractions beyond MS whereas this is not necessary for the diagrams (7). The subtraction 
term for the Dirac form factor reads 



5 In order to take care of the mass shift due to C^jj, we use to 2 — m — Ac\M 2 in the fermion propagators 
At the end of the calculation we then replace ni2 by the physical nucleon mass rnjv, because the 
difference shows up only in higher orders. 
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and similarly for the Pauli form factor 6 



2 

AP 5 9°A m N (m - 4ciM 2 ) 

AF 2 = (3-7S), 

2 

8 g° A m N (m - 4ciM 2 ) 
= 8^ T3 ' 

AF 9 » = - ^^^f 1 ^ (3c 7 -2c 6 r 3 ). 

For completeness, we determine the expansion of the couplings of the MS scheme Cg, c 7 , dg, 
d?, eg 4 , and e 74 in terms of our renormalized parameters, 7 



o 2 

9a rn 
64tt 2 F 2 

2 

3 g A m 



it A no . 



c 7 = C7 + 327 r 2 F 2 ( 1 + 2mc? )- 

o 2 

at a 9a m 
~ de ~ 32^ C6 ' 

o ° 2 

^ = ^ + 128^ C7 ' 



Q ° 2 

3 



e 74 = 674 + 7TZ^2 C 6' 



654 " 654 2 5 67r 2 F 2C7 ' 

o 2 

64rr 2 F 2 ' 

o 2 

e^os = e i05 + 1 ?^ 2 1 (3 + 12mc 7 ), 



'106 



128tt 2 F 2 

o 2 
fi'A Ci 

e i06-^^(5 + 8mc 6 ), 



producing all the required counterterm contributions. 

In order to obtain the final expression one has to sum up all contributions and to multiply 
the result with the wave function renormalization constant Z, 8 given by 



9g A M 
32tt 2 F 2 



2 o o 2 .> 

" ' 2 r 2 /M x 1 ° " ■ n/r3 



+ In , 

3 \m 



+ ^M1 + (M% (11) 
647rF 2 m v y v 7 



Next, we need to fix the parameters in order to get a graphical representation of the form 
factors. The values of the renormalized LECs c 2 and c 4 are taken from a (tree-level) fit [37| 



6 One factor of tun is due to the fact that the parameterization of Eq. JBJ contains an explicit factor 1/rnjv 
multiplying the Pauli form factor. 

7 For notational convenience, we omit a subscript R for the EOMS-renormalized parameters, but one should 
not confuse them with the bare parameters of the bare Lagrangian. 

8 For F 2 , one needs Z only up to order 0{q 2 ) due to the factor q v in Eq. ©. 
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to the irN scattering threshold parameters of Ref. 38]: C2 = 2.5m^ 1 and C4 = 2.3m^ r 1 . Since 
the q only enter loop expressions at 0(q 4 ) a determination at 0(q 2 ) is consistent with the 
present accuracy. 

The chiral expansion of the anomalous isoscalar and isovector magnetic moments can be 
written as 



(s/v) 



(s/v) 



+ k[ s/v) M 



.(«/«) 



/ >2 ' .I/" 111 ( J + /,';, 



(s/v) 



M 2 + 0(M 3 ), 



(12) 



where the coefficients, in terms of the EOMS-renormalized parameters, are given by 



K 



.(«) 
'0 



(«) 



2 07m, 
0, 

3(1 + 2c 7 m) g° A 
8tt 2 F 2 

o 2 

3(2 + 3c 7 m) qa 



16tt 2 F 2 



— 8ciC 7 — 32eio5?7t, 



(13) 



k q 

» 



(v) 
2 



o 2 

9 a rn 
4ttF 2 ' 



o 2 



4(c 4 - c 6 )m- g>A (7 + 8c 6 m) 
8vr 2 F 2 
/a (l + 5c 6 m) 
8vr 2 F 2 



16ciC 6 — 16eio6^- 



(14) 



Comparing with Ref. [33j we see that the lowest-order terms as well as the non-analytic 
terms ~ M ~ m 1 / 2 and ~ ln(M) coincide, while the analytic K3 terms differ, because we 
use a different renormalization scheme. 

Since chiral perturbation theory does not predict the anomalous magnetic moments, we 
take the empirical values as obtained from Eq. (JJJ) to fix the Pauli form factors at Q 2 = 0. 
The constants d 6 , dj, e^, and eu were obtained from the charge and magnetic radii of the 
nucleon, which are defined via the Sachs form factors [see Eqs. Q and ([10))]: 



((r p E ) 2 ) 



{{r n E ?) 



6 dG p E (Q 2 



G P E (0) dQ 2 
dGUQ 2 ) 



Q 2 =0 



-6 



dQ 2 



((r n M f) 



6 dG p M (Q 2 ) 



G P M (0) dQ 2 
6 dG n M (Q 2 



Q 2 =0 



Q 2 =0 



G n M (0) dQ 2 



For the numerical analysis, we make use of the charge and magnetic radii as obtained in the 



and = 0.889 fm. Within the accuracy of our calculation the constants cq, C7, eios, and 



dispersion-theoretical analysis of Ref. [2 
and r n M = 0.889 fm. Within the 
eio6 appear in the combinations 

c 6 = c 6 -4M 2 ei 06 , 



0.847 fm, t% 



-0.113 fm, r p M 



0.836 fm 



c 7 = c 7 - l6M 2 e W5 . 



(15) 



Therefore for the numerical analysis we only need the values for the combinations cq^. The 
values for the LECs are summarized in Table UJ From Table H] we see that there is quite a big 
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difference in the value for the constant de as obtained in the infrared regularization and our 
renormalization scheme, respectively. The constants d§ and d-j yield the leading order of the 
electromagnetic radii and thus determine, together with those loop graphs which depend on 
Q 2 , the slope of the graphs at the origin. 





C2 


c 4 


C6 


C7 


d 6 


dj 


e54 


e 74 


EOMS 


2.66 


2.45 


1.26 


-0.13 


-0.69 


-0.50 


0.19 


1.59 


IR 


2.66 


2.45 


1.26 


-0.18 


0.54 


-0.73 


0.25 


1.93 



TABLE I: Values for the relevant low-energy constants in the extended on-mass-shell renormal- 
ization scheme (EOMS) and in the infrared regularization scheme (IR) within an 0(q 4 ) calculation 
of relativistic chiral perturbation theory. The LECs q are given in units of GeV -1 , the di in units 
of GeV -2 , and the ei in units of GeV -3 . 

The results for the Dirac and Pauli form factors are shown in Fig. El for momentum 
transfers < Q 2 < 0.4 GeV 2 . For comparison the results of the infrared regularization [33] 
are also shown. We see from Fig. El that for the Dirac form factor of the proton and for 
both Pauli form factors the two methods generate very similar results. However, this is not 
the case for the Dirac form factor of the neutron which is very small in the entire small Q 2 
region. We conclude that F™ significantly depends on higher-order contributions. 

As an example, let us have a closer look at the Dirac form factor of the proton, Ff , which, 
for both methods, shows an almost linear behavior in Q 2 , i.e., has very little curvature. In 
order to show the influence of the loop contributions, in Fig. |3] we compare the full results 
for Ff in both relativistic renormalization schemes with the contact-graph contributions 
only. We see that the loop contributions cannot be neglected, i.e., the contribution of 
the pion cloud plays an important role. Moreover, in the infrared regularization the main 
contribution to the slope is due to loop contributions, while in our renormalization scheme 
the loop contributions are smaller. The situation here is different from the electromagnetic 
form factor of the pion whose linear behavior is predominantly due to a contact graph 0. 



C. Results for the Sachs form factors 

We now turn to a discussion of the Sachs form factors in the EOMS renormalization 
scheme and compare them with experimental data as well as the results of the infrared 
regularization and of HBxPT at 0(q 3 ) (HUD!- 9 Recall that a full 0(q A ) calculation yields 
Ge and Gm up to 0(q 3 ) and 0(q 2 ), respectively. 10 

In order to discuss the Sachs form factors resulting from an 0(q 3 ) calculation, we need 
to recalculate the adapted values for the LECs (see Table HTJ) since the values of Table 
were obtained from an 0(q 4 ) calculation. Note that the changes in the LECs turn out 
to be reasonably small. Figure H] shows our results for the Sachs form factors at 0(q 3 ) 
together with experimental data as well as the infrared regularization and HB%PT results 



9 The effects due to t he A resonance, as caicuiated in the smaii-scaie expansion, can be found in Ref. 32]. 
As discussed in Ref. |33J, when combining Fi and F 2 to give the Sachs form factors Ge and Gm [see Eqs. 
10 and (|10l) ] the chiral orders are mixed. 
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C6 


C7 


da 


d-j 


EOMS 
IR 


1.34 
1.36 


-0.15 
-0.23 


-0.70 
0.50 


-0.49 
-0.73 



TABLE II: Values for the relevant low-energy constants in the extended on-mass-shell renormal- 
ization scheme (EOMS) and in the infrared regularization scheme (IR) within an 0(q 3 ) calculation 
of relativistic chiral perturbation theory. The LECs Cj are given in units of GeV -1 and the di in 
units of GeV -2 . 

at 0(q 3 ). For both Sachs form factors of the proton and the magnetic Sachs form factor of 
the neutron the two relativistic calculations are almost identical while the HB^PT results 
clearly deviate from the relativistic calculations. This difference may be interpreted as being 
due to relativistic corrections which, in HB^PT, would have to show up in higher orders. It is 
somewhat surprising that in all three cases the HB^PT results are closer to the experimental 
data. The slopes of the electric form factors are essentially determined by fitting the mean 
square electric radii. On the other hand, at 0(q 3 ) the magnetic radii do not contain any 
free parameters. This explains the difference in the behavior of the slopes of the magnetic 
form factors between the relativistic calculation on the one hand and the HB%PT calculation 
on the other hand, originating from the fact that the relativistic calculations at 0(q 3 ) for 
the magnetic radii also contain higher-order terms in M, whereas HB%PT only involves the 
leading-order term. Finally, the electric form factor of the neutron is very small and all three 
theoretical curves disagree. Here, the relativistic 0(q 3 ) calculations are in better agreement 
with the experimental data. 

Figure El shows the results of the two relativistic calculations for the Sachs form factors at 
0(q A ). The description of G P E , G P M , and G\j is only marginally better than that of the 0(q 3 ) 
calculation. For the very-small Q 2 region this is due to additional free parameters which have 
been adjusted to the magnetic radii. Clearly, the 0(q 4 ) results do not generate sufficient 
curvature which is emphasized in Figs. El El and|Hl showing the form factors and data divided 
by the dipole form factor Go and normalized to one at the origin. Finally, the situation for 
G E seems to be even worse, where we find better agreement with the experimental data for 
the 0(q 3 ) results. We conclude that the perturbation series converges, at best, slowly and 
that higher-order contributions must play an important role. 

IV. SUMMARY 

We have calculated the electromagnetic form factors of the nucleon at 0(q A ) (one- loop 
order) in relativistic chiral perturbation theory using the extended on-mass-shell renormal- 
ization scheme of Ref . [I?) • The relevant low-energy coupling constants have been fitted to 
the anomalous magnetic moments and to the charge and magnetic mean square radii. As 
a general trend the results for both the proton form factors G P E and G P M and the magnetic 
neutron form factor G 7 ^ do not show sufficient curvature to generate agreement with the 
experimental data beyond small values of Q 2 ~ 0.1 GeV 2 . Our results are very similar to 
those of the infrared regularization [3j| and the small differences between the two methods 
are related to the way how the regular higher-order terms of loop integrals are treated. In 
the case of the electric neutron form factor G E , in both renormalization schemes the de- 
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scription does not seem to improve when going from 0(q 3 ) to 0(q 4 ). We conclude that a 
relativistic treatment at the one-loop level using nucleon and pion degrees of freedom only 
is not sufficient to describe the form factors for Q 2 > 0.1 GeV 2 and that higher-order con- 
tributions must play an important role. Such higher-order contributions either have to be 
studied in the framework of a two-loop calculation or by explicitly including other dynamical 
degrees of freedom such as vector mesons [33|. We stress that the EOMS renormalization 
scheme allows for also implementing a consistent power counting in relativistic baryon chiral 
perturbation theory when vector (and axial- vector) mesons are explicitly included [56j]- 
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APPENDIX A: LOOP INTEGRALS 

1. Definition of the loop integrals 

In the appendix we use the following notation 

P»=tf+lf f , q" = P>}-rf, t = q 2 . 
The loop integrals with one and two internal lines are defined as follows: 

. , d n k 1 



I N = iJ 



g*I$(t)=i 



tgriW(t)+<f<fl!*p(t)=i 

lNN(t) = i 

I-kn{v 2 ) = i J 



(2tt)" k 2 - M 2 + it 

d n k 1 
(27r) n fc 2 -m 2 + ie 
d n k 1 
(2tt)™ [k 2 -M 2 + ie][{k + q) 2 — M 2 + ie] ' 
d n k k^ 
(2n) n [k 2 — M 2 + ie][(k + q) 2 - M 2 + ie] 
d n k k»k v 



(2tt)" [k 2 — M 2 + ie] [(k + q) 2 — M 2 + ie] 
d n k 1 
(27r) n [k 2 -m 2 + ie] [(k + q) 2 - m 2 + ie}' 
d n k 1 
(2ir) n [k 2 — M 2 + ie] \{k - p) 2 - m 2 + ie] ' 
d n k k» 



(2tt)" [k 2 — M 2 + ie] [(k - p) 2 - m 2 + ie}' 

In the following loop integrals involving three internal lines we always assume on-shell kine- 
matics, pi = p] = m 2 N) 

f d k 

lmTN{t) = I 



(2tt)" [k 2 -M 2 + ie] [(k + q) 2 — M 2 + ie] [{k - Pl ) 2 - m 2 + ie] 
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P*l£> N {t)-±q*I mIf {t)=iJ 



d n k 



k 4 ' 



(2tt)™ [k 2 — M 2 + ie] [(fc + q) 2 — M 2 + z'e] [(jfe - p*) 2 - m 2 + ie] 



d n k 



k»k u 



InNN(t) — i 



(2n) n [k 2 — M 2 + ie] [(k + q) 2 — M 2 + ie] [(jfe - p { ) 2 - m 2 + ie] ' 
d n k 1 

(2tt)™ [A; 2 — M 2 + ie] [{k - Pi ) 2 - m 2 + ie] [(k - p f ) 2 - m 2 + ie] ' 



d n k 



d n k k^k u 



(2yr) n [k 2 — M 2 + ie] [(Jfe - p^ 2 - m 2 + ie] [(jfe - p/) 2 - m 2 + ie] ' 



r(W), 



(27r) n [A; 2 — M 2 + ie] [(jfe - p t ) 2 - m 2 + ie] [(jfe - p/) 2 - m 2 + ie] ' 



2. Scalar loop integrals 

Defining 



A 



m 



n— 4 



1 1 



and 



16tt 2 U-4 2 

p 2 -m 2 - M 2 



-[ln(4vr) + r'(l) + l] , 



n 



the scalar loop integrals are given by (l7| 



2mM 



where 



„- M 2 /M\ 
J^ = 2M 2 A + — In — , 
07r z Km J 



in 



2m 2 A, 



I nn (t) = 2 A + 



/^(p 2 ) =2A + 



167T 2 

/jviv(t) = 2A + 
1 



I ,2,n|-| +J (»)( ]j L) 



m / 



1 



t 



1 + J (0) f — 



16tt 2 



-1 



167T 2 

p 2 -m 2 + M 2 (M\ 2mM 



In 



+ 



P 



F(Q) 



J (0 \x) = f dz\n\l + x(z 2 - z) -ie] 

JO 



x < 0, 



-2 + 2^1 - 1 arccot (y^ - l) , < x < 4, 
-2 -crln(i=^) -i7rcr, 4 < x, 
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with 



and 



a(x) = Jl--, x£[0,4], 



Vft 2 - 1 In [-n - Vft 2 - 1 



F{Q) = { Vl-ft 2 arccos(-ft) 

Vft 2 - 1 In (fi + Vft 2 - 1 



ft < -1, 
-1< n < 1, 



?7T 



Vft 2 - l, l < ft. 



In our numerical analysis, we made use of a parameterization in terms of a Feynman pa- 
rameter integral to analyze the loop integrals with three internal lines. Such integrals may 
also be determined using a dispersion relation (see App. C of Ref. 

3. Reduction of the tensorial loop integrals 

The reduction of the tensorial loop integrals to the corresponding scalar ones can be 
performed in the standard way |57[. We obtain the following results: 

4S(*) = 



(n - l)t 
1 

7 



n — 2 T (n , , 9 \ r . . 

-4 + -(t- M 2 ) I n Jt) + — ( 3M 2 - - 
3 3 v / mW 144vr 2 V 2 



4° 0) (*) 



(n - l)t 
1 

7 6~" w ' 12 



I/ 7r + I(4M 2 -t)/ 7r7r (t) 



^ 7r (t) + -^(4M 2 -t)/_(t) 



1 / M 2 t 



In 2 V 6 36. 



2p 2 



2(4m 2 N - t) 



I N -h + (p 2 -m 2 + M 2 )^ N (p 2 
(2M 2 - t)I^ N {t) + 2/^(7^) - 2/ OT (t) 



r(00) (f] 



-^—^ [2h N {0 + (4M 2 - t)I nwN (t) - 2(2M 2 - t)li P J N (t) 
I \2h N {m 2 N ) + (4M 2 - t% wN (t) ~ 2(2M 2 - t)I { n P J N (t) - -L 

o L o7T . 



7 (PP) M 



4(n-2)(4ra 2 v - t) 



-2/ wJV (m^) + 2(n - 2)I^(m 2 N ) - (4M 2 - t)I n7TN (t) 
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+2(n-l)(2M 2 -t)I^ N (t) 
1 



8(4m% - t) 



2I nN (m 2 N ) + 4lH(m 2 N ) - {AM 2 - t)I wnN (t) 



+6(2M 2 -t)lZUt) + 



8tt 2 



T (<?<?) 
1 ttttN 



(t) 



tN^n) 



2(n - 2)J$ «) 



_ 4M 2 -{n- l)t] hMt) + 2(2M 2 - t)I { /J N (t)} 



1 

8t 



2/ ffJV (m^) - AI^{m 2 N ) - (4M 2 - 3t) Wt) + 2(2M 2 - t)I { ^ N {t) 



+ 



T {P) 



(t) 



1 



4m 2 - t 



r (oo) 

ttNN 



1 



ra-2 
1 



M 2 I irNN (t) - I nN (m 2 N ) + I NN (t) 



WW - M 2 + ^/ JViV (t) 



WW-/SjvW M 2 + ^/^(t) - 



32tt 2 J ' 



8tt 2 



ttNN 



(*) 



(n-2)(Am 2 N -t) 



n - l)lfJ N (t) - I w NN{t) 



, r2 n—2 T 
M Z 



(p) / 2 \ i n 3 
ttNv^n) H 7^— *nn 



( 



2(4m 2 r -0 



(P) ^2 



32tt 2 



/■(ot) 

ttNN 



(*) 



(n - 2)t 



M 2 + - ± WO 



2t 



- W(t)J M 2 + - - W*) + ^ j 



(P) f™ 2 



1 1 



APPENDIX B: CONTRIBUTIONS TO THE DIRAC AND PAULI FORM FAC- 
TORS 

The unrenormalized expressions for the Dirac form factor F\ read 

,1+3 + n 



- (r 3 d 6 + 2d 7 ) t, 



F? = 



;3 - r 3 ) {I n + Am 2 N I^{m 2 N ) - Am 2 N [M 2 I nNN (t) + I NN (t)} 



o 2 
9 A 

8F 2 

+8m 2 N I%(t) + 32m%li P N I 2(t)}, 

,6 = T 3 

2F 2 
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2 

F 1 7 = ^T 3 [I« + 2mllM(mU- 



o 2 



The contributions to the Pauli form factor F 2 are given by 

F 2 2+3+4 = 2m N T 3 c 6 + m^cy + (r 3 cf 6 + 2cf 7 ) £ + 2m N (2 e 54 + r 3 e 74 ) t 
-8m N M 2 (2 eio5 + r 3 ei 06 ) , 

^ 5 = -^(3-r 3 )4m^/S(t), 

o 2 



F 2 8 = ^r 3 16myS(t), 



o 2 



p io = (3C7 _ 2c6r3 ) {/ T + H/SK) + HM 2 

+Am%I NN (t) - IQrnlCUt) + 8m% t [^(t) - 
ii 3M 2 , / M 2 \ m N 

The upper indices refer to the corresponding diagrams of Fig. ^ Note that we did not 
distinguish between m and mjv, because the difference only shows up at higher orders. Our 
results agree with the expressions of Ref. up to terms vanishing in infrared regularization. 



[1] S. Weinberg, Physica A 96, 327 (1979). 

[2] J. Gasser and H. Leutwyler, Annals Phys. 158, 142 (1984). 

[3] S. Scherer, in Advances in Nuclear Physics, Vol. 27, edited by J. W. Negele and E. W. Vogt 

(Kluwer Academic/Plenum Publishers, New York, 2003). 

[4] J. Gasser, M. E. Sainio, and A. Svarc, Nucl. Phys. B307, 779 (1988). 

[5] E. Jenkins and A. V. Manohar, Phys. Lett. B 255, 558 (1991); ibid. 259, 353 (1991). 

[6] V. Bernard, N. Kaiser, J. Kambor, and U.-G. Meifiner, Nucl. Phys. B388, 315 (1992). 

[7] V. Bernard, N. Kaiser, and U.-G. Meifiner, Nucl. Phys. A611, 429 (1996). 



15 



[8] G. Ecker and M. Mojzis, Phys. Lett. B 365, 312 (1996). 

[9] N. Fettes, U.-G. Meifiner, M. Mojzis, and S. Steininger, Annals Phys. 283, 273 (2000); ibid. 
288, 249 (2001). 



[10] H. Tang, hep-ph/9607436 



[11] P. J. Ellis and H. B. Tang, Phys. Rev. C 57, 3356 (1998). 
[12] T. Becher and H. Leutwyler, Eur. Phys. J. C 9, 643 (1999). 
[13] J. Gegelia and G. Japaridze, Phys. Rev. D 60, 114038 (1999). 
[14] J. Gegelia, G. Japaridze, and X. Q. Wang, |hep^ph/9910260| 
[15] M. Lutz, Nucl. Phys. A677, 241 (2000). 

[16] M. F. Lutz and E. E. Kolomeitsev, Nucl. Phys. A700, 193 (2002). 
[17] T. Fuchs, J. Gegelia, G. Japaridze, and S. Scherer, hep-ph/0302117, 
[18] F. J. Ernst, R. G. Sachs, and K. C. Wali, Phys. Rev. 119, 1105 (1960). 
[19] R. G. Sachs, Phys. Rev. 126, 2256 (1962). 
[20] J. J. Kelly, Phys. Rev. C 66, 065203 (2002). 
[21] J. Friedrich and Th. Walcher, lhep^ph/0303054. 



[22] F. Iachello, A. D. Jackson, and A. Lande, Phys. Lett. B 43, 191 (1973). 

[23] P. Kroll, M. Schurmann, and W. Schweiger, Z. Phys. A 338, 339 (1991). 

[24] X. D. Ji, Phys. Lett. B 254, 456 (1991). 

[25] F. Schlumpf, J. Phys. G 20, 237 (1994). 

[26] F. Cardarelli, E. Pace, G. Salme, and S. Simula, Phys. Lett. B 357, 267 (1995). 

[27] C. V. Christov, A. Z. Gorski, K. Goeke, and P. V. Pobylitsa, Nucl. Phys. A592, 513 (1995). 

[28] P. Mergell, U.-G. Meifiner, and D. Drechsel, Nucl. Phys. A596, 367 (1996). 

[29] D. H. Lu, A. W. Thomas, and A. G. Williams, Phys. Rev. C 57, 2628 (1998). 

[30] S. Both, L. Y. Glozman, W. Klink, W. Plessas, M. Radici, and R. F. Wagenbrunn, Eur. Phys. 
J. A 14, 17 (2002). 

[31] H. W. Fearing, R. Lewis, N. Mobed, and S. Scherer, Phys. Rev. D 56, 1783 (1997). 

[32] V. Bernard, H. W. Fearing, T. R. Hemmert, and U.-G. Meifiner, Nucl. Phys. A635, 121 

(1998); ibid. A642, 563 (1998). 

[33] B. Kubis and U.-G. Meifiner, Nucl. Phys. A679, 698 (2001). 

[34] T. Ebertshauser, H. W. Fearing, and S. Scherer, Phys. Rev. D 65, 054033 (2002). 

[35] G. Colangelo, J. Gasser, and H. Leutwyler, Phys. Rev. Lett. 86, 5008 (2001). 

[36] S. Scherer, A. Y. Korchin, and J. H. Koch, Phys. Rev. C 54, 904 (1996). 

[37] T. Becher and H. Leutwyler, JHEP 0106, 017 (2001). 

[38] R. Koch, Nucl. Phys. A448, 707 (1986). 

[39] L. E. Price, J. R. Dunning, M. Goitein, K. Hanson, T. Kirk, and R. Wilson, Phys. Rev. D 4, 
45 (1971). 

[40] C. Berger, V. Burkert, G. Knop, B. Langenbeck, and K. Rith, Phys. Lett. B 35, 87 (1971). 

[41] K. M. Hanson, J. R. Dunning, M. Goitein, T. Kirk, L. E. Price, and R. Wilson, Phys. Rev. D 
8, 753 (1973). 

[42] G. G. Simon, C. Schmitt, F. Borkowski, and V. H. Walther, Nucl. Phys. A333, 381 (1980). 

[43] T. Eden et al, Phys. Rev. C 50, 1749 (1994). 

[44] I. Passchier et al, Phys. Rev. Lett. 82, 4988 (1999). 

[45] M. Ostrick et al, Phys. Rev. Lett. 83, 276 (1999). 

[46] C. Herberg et al, Eur. Phys. J. A 5, 131 (1999). 

[47] J. Becker et al, Eur. Phys. J. A 6, 329 (1999). 

[48] T. Janssens, R. Hofstadter, E. B. Hughes, and M. R. Yearian, Phys. Rev. 142, 922 (1966). 



16 



[49] G. Holder, E. Pietarinen, I. Sabba Stefanescu, F. Borkowski, G. G. Simon, V. H. Walther, 

and R. D. Wendling, Nucl. Phys. B114, 505 (1976). 
[50] P. Markowitz et al, Phys. Rev. C 48, 5 (1993). 
[51] H. Anklin et al, Phys. Lett. B 336, 313 (1994). 
[52] E. E. Bruins et al, Phys. Rev. Lett. 75, 21 (1995). 
[53] H. Anklin et al, Phys. Lett. B 428, 248 (1998). 
[54] W. Xu et al, Phys. Rev. Lett. 85, 2900 (2000). 
[55] G. Kubon et al, Phys. Lett. B 524, 26 (2002). 
[56] T. Fuchs, J. Gegelia, and S. Scherer, in preparation. 
[57] G. Passarino and M. J. Veltman, Nucl. Phys. B160, 151 (1979). 



17 



(1) 



(2) 



(3) 



(4) 





(10) (11) (12) 



FIG. 1: Feynman diagrams contributing to the electromagnetic form factors up to and including 
0(q A ). External-leg corrections are not shown. Solid, dashed, and wiggly lines refer to nucleons, 
pions, and photons, respectively. 
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FIG. 2: The Dirac and Pauli form factors of the nucleoli in relativistic chiral perturbation theory 
at 0(<? 4 ). Full lines: our results in the extended on-mass shell (EOMS) scheme; dashed lines: 
infrared regularization result |33| . 
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Q 2 [GeV 2 ] 



FIG. 3: The Dirac form factor of the proton at 0(q 4 ). Solid line: our result; dashed line: infrared 
regularization; dotted line: our result without loop contribution; dashed-dotted line: infrared- 
regularization result without loop contribution. 
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FIG. 4: The Sachs form factors of the nucleoli at 0(q 3 ). The solid, dashed, and dotted lines 
refer to the results in the EOMS scheme, the infrared regularization ji^], and HByPT 



resp ectively. The experimental data for G P E: Crj L , a nd G™, are taken from Refs. 



resp ectively, ine experimental data lor &' F , U b , (j ^, a 

20,0,0,143, 0,00,0, and OEJEaEl 




31], 



5J, |55|, respectively. 
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FIG. 5: The Sachs form factors of the nucleon at 0(q 4 ). The solid and dashed lines refer to the 
results in the EOMS scheme and the infrared regularization 33 1, respecti vely . The experimental 



data for G|, GPl, G p t 



and CI, are taken from Refs. 



3 , and 0E], S El, EJ,I3 , respectively. 



zation 6o\, respectively, ine experime 

30, HQ, BHHHQ, 
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FIG. 6: The magnetic form factor of the neutron divided by The solid and dashed lines 

refer to the results in the EOMS scheme at 0(q 4 ) and 0((/ 3 ), respectively. The dotted line is the 
linear approximation of the dipole form factor, i.e., [1 — 2Q 2 /(0.71GeV 2 )]/Gz)((5 2 )- The data are 

H3, E3, E3, E3, E3, E3i 



taken from Refs. 
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FIG. 7: The magnetic form factor of the proton divided by \x v Gr>. The solid and dashed lines 
refer to the results in the EOMS scheme at 0(q 4 ) and 0((/ 3 ), respectively. The dotted line is the 
linear approximation of the dipole form factor, i.e., [1 — 2Q 2 /(0.71GeV 2 )]/G/)((5 2 ). The data are 

taken from Refs. 0,0, SB. 
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FIG. 8: The electric form factor of the proton divided by Go- The solid and dashed lines refer 
to the results in the EOMS scheme at 0(q 4 ) and 0(q i ), respectively. The dotted line is the linear 
approximation of the dipole form factor, i.e., [1 — 2Q 2 /(0.71GeV 2 )]/Gn(Q 2 )- The data are taken 
from Refs. 



25 



